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Single molecule magnets and single spin centers can be individually addressed when coupled 
to contacts forming an electrical junction. In order to control and engineer the magnetism of 
quantum devices, it is necessary to quantify how the structural and chemical environment of the 
junction affects the spin center Metrics such as coordination number or symmetry provide 

a simple method to quantify the local environment, but neglect the many-body interactions of an 
impurity spin when coupled to contacts [7]. Here, we utilize a highly corrugated hexagonal boron 
nitride (h- BN) monolayer 8]0J to mediate the coupling between a cobalt spin in CoH x (x = 1,2) 
complexes and the metal contact. While the hydrogen atoms control the total effective spin, the 
corrugation is found to smoothly tune the Kondo exchange interaction between the spin and the 
underlying metal. Using scanning tunneling microscopy and spectroscopy together with numerical 
simulations, we quantitatively demonstrate how the Kondo exchange interaction mimics chemical 
tailoring and changes the magnetic anisotropy. 

*P. J. and T. H. contributed equally. 

Magnetic anisotropy defines the stability of a spin in a 
preferred direction [2]. For adatoms on surfaces, the low 
coordination number and changes in hybridization can 
lead to dramatic enhancement of magnetic anisotropy 
[3 Ho]. Surface adsorption site and the presence of hy¬ 
drogen has been shown to alter the magnetic anisotropy 
of adatoms on bare and graphene covered Pt(lll) fill 
H3j . Furthermore, the exchange interaction and strain 
has been invoked for 3d adatoms on CU 2 N islands where 
the adatom position on the island affects the observed 
magnetic anisotropy mm- Studies on single molecule 
magnets (SMMs) containing 3d or 4/ spin centers have 
revealed that chemical changes to the ligands surround¬ 
ing the spin affect the magnetic anisotropy m- How¬ 
ever, the most important factor for maintaining magnetic 
anisotropy in SMMs is a low coordination number and a 
high axial symmetry HH3CLZ!. 

Magnetic anisotropy is not guaranteed in SMMs or sin¬ 
gle spin centers upon coupling to contacts. The spin in¬ 
teracts with the electron bath through the exchange in¬ 
teraction leading to a finite state lifetime and the decay 
of quantum coherence PEQ!. Additionally, the scatter¬ 
ing of the spin with the electron bath results in an en¬ 
ergy renormalization of the spin’s eigenstate energy lev¬ 
els, similar to the case of a damped harmonic oscillator 
[18 . In practice, this leads to a net reduction of the mag¬ 
netic anisotropy, pushing the system closer to a Kondo 
state. At the heart of the Kondo effect are spin-flip scat¬ 
tering processes between localized states at the impurity 
spin and delocalized states in the bulk conduction band, 
resulting in the formation of a correlated quantum state 
[20] . The Kondo regime is reached when the magnetic 
moment of the impurity spin is screened by the electron 


bath, with the exchange interaction defining the relevant 
energy scale, the Kondo temperature (TV). High spin 
systems with a total spin S > 1/2 have the potential for 
both magnetic anisotropy and the Kondo effect [21] 22] . 
Thus, the Kondo exchange interaction with the electron 
bath can force the impurity spin into a competing Kondo 
state, where antiferromagnetic coupling with the reser¬ 
voir reduces or even quenches the magnetic moment. The 
outcome of this competition can be determined in local 
transport measurements, but few quantitative measures 
of this competition exist. 

Here, we study CoH x complexes coupled to a spa¬ 
tially varying template, the /i-BN/Rh(lll) moire, to ob¬ 
serve and model how the environment influences mag¬ 
netic anisotropy. The h -BN monolayer, a wide band gap 
two dimensional material, decouples and mediates the 
interactions between CoH x and the underlying Rh metal 
while lattice mismatch leads to a spatial corrugation re¬ 
sulting in an enlarged unit cell with 3.2 nm periodicity 
corresponding to 13 BN units on top of 12 Rh atoms 
[8]. The local adsorption configuration of CoH x on the 
h -BN is conserved across the moire unit cell, with the 
large number of inequivalent adsorption sites allowing us 
to explore how hybridization affects magnetic anisotropy. 
To complement our experimental observations, we model 
transport through the CoH x complexes using Hamilto¬ 
nians that incorporate magnetic anisotropy as well as 
coupling to the environment. This is accomplished by 
parameterizing the environment through use of a dimen¬ 
sionless coupling constant — Jpo, describing the strength 
of the Kondo exchange interaction, J, between the lo¬ 
calized spin and the electron density po of the substrate 
near the Fermi level (see SI). 
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Figure 1. CoH^ adsorbed on a h— BN/Rh(lll) surface. 

(a) Constant current STM topography with three CoB x com¬ 
plexes (protrusions) adsorbed on different sites (15 x 7 nm 2 
image size, V = —100 mV, I — 20 pA, T = 1.4 K). High sym¬ 
metry points of the moire are marked by the white overlay. 

(b) Sketch of the atom positions for the adsorption of CoH. 
The h -BN registry with Rh(lll) shifts across the moire unit 
cell with three high symmetry sites: at the valley site (v) the 
Rh is directly underneath the N, whereas for the two unequal 
rim sites (n and r 2 ) changes in the registry and distance to 
the surface are observed, (c) Line profiles along the dashed 
lines indicated in (a) show two CoH x systems with adsorption 
sites ri and v (red line) and a /i-BN reference cut (blue line, 
offset by 0.5 A), (d) Differential conductance dl/dV curves 
versus bias voltage of three different CoH x systems (stabi¬ 
lization setpoint: I = 500 pA, V = —15 mV, T = 1.4 K, 
curves vertically offset for clarity). The upper curve (grey) 
shows a spin-1/2 Kondo resonance centered at zero bias. The 
two lower curves (red and blue) show step-like conductance 
increases symmetric around zero bias indicating a spin-1 sys¬ 
tem. Solid black lines are least-square fits using a perturbative 
transport model. 


Figure shows a representative scanning tunnel¬ 
ing microscopy (STM) topograph of the /i-BN/Rh(lll) 
moire with isolated CoH x (x = 1,2) complexes, line pro¬ 
files across the h -BN indicate CoH x can adsorb at mul¬ 
tiple positions within the moire (Figure E>) m- On 
these CoH x complexes we measure the differential con¬ 
ductance, dl/dV , against the applied bias voltage V be¬ 
tween tip and sample at low-temperature (T = 1.4 K) 
and zero magnetic field {B = 0 T, details see Methods). 
The spectra can be divided into two broad classes: a 
sharp peak centered at zero bias or two symmetric steps 
of increasing conductance at well-defined threshold en¬ 
ergies (Figure [lj:). The peak at zero bias is consistent 
with a spin-1/2 Kondo resonance while the steps corre¬ 
spond to the onset of inelastic excitations from the mag¬ 
netic ground state to excited states. The observation of 
two steps hints at a spin-1 system with zero field split¬ 
ting. The two lower spectra (Figure [1J; red, blue curves) 
are measured on CoH at different parts of the moire and 
share the same characteristics but the step positions vary. 

We employ density functional theory (DFT) to corre¬ 
late the magnetic properties of CoH x with the local ad¬ 
sorption configuration. Our calculations (see Methods) 
show that adsorption in the BN hexagon, i.e. hollow 
site, is preferable for bare Co. The addition of hydrogen 
shifts the preferred adsorption site to N, with the hollow 
site adsorption energy consistently higher. For CoH com¬ 
plexes the preferred hydrogen position was found to be 
either exactly on top of Co or tilted towards the nearest 
B atom (Figure |2^i). An important consequence of the 
N adsorption site is the linear crystal field acting on the 
cobalt {i.e. N — Co — H) removing the 5-fold degeneracy 
of the d -levels (Figure^). 

In Figure^ the spin-resolved, symmetry decomposed 
local density of states of CoH and C 0 H 2 adsorbed in the 
h -BN valley is plotted. The atomic d-levels are split 
roughly 1.2 eV by the intrinsic Stoner exchange giving a 
bare Co adatom a magnetic moment of 2.2 Bohr magne¬ 
tons (/jLb)- Formation of CoH leads to hybridization of 
the H sp orbitals and the Co orbitals, slightly reducing 
the magnetic moment to 2.0 /ip, equivalent to a 3d 8 con¬ 
figuration (Figure^). The second hydrogen changes the 
picture significantly, with the sp — d hybridization suffi¬ 
cient to bring the Co d-levels closer together, reducing 
the magnetic moment to 1.2 /ip resulting in a 3d 9 config¬ 
uration. Therefore, from our spectroscopic observations 
and DFT calculations we identify CoH as an effective 
spin-1 and C 0 H 2 as spin-1/2 system. 

Figure [2]i shows the spin density distribution for CoH 
in a N adsorption configuration at two high symmetry 
points of the moire. The strong vertical bond between 
Co and N leads to an effective spin-polarization along 
this axis and can be expected to provide the system with 
out-of-plane magnetic anisotropy. Tilting of the hydro¬ 
gen and the underlying lattice mismatch reduces the Cs v 
symmetry and introduces small shifts in the d xz , d yz lev- 







3 


a 



b 


Z2 


A l 


xz.yz+Jfr- I _ 

xyX-y’ -H- “H“ \ 2 [E 


10) 

ID+l-l) 
|l>—l-l > 



d 



|-0.5 Polarizatioi7(%r^^^B 


O Co ° H • B ON ® Rh 


Figure 2. CoH and C0H2 density of states, (a) Ball and stick model of the adsorption of CoH on h- BN. The linear 
adsorption geometry of CoH on the N atom is emphasized and marks the main (axial) magnetic anisotropy (D) along the 
z-axis. Additional transverse anisotropy (E) in the x — y plane further breaks the symmetry, (b) Schematic linear crystal field 
splitting diagram for the 3 d 8 shell of Co highlighting the origin of the axial (D) and transverse ( E ) magnetic anisotropy. The 
magnetic ground state is an antisymmetric superposition of mn z — |+1) and | —1) states ( mn z is the magnetic moment in units 
of the reduced Planck constant h in ^-direction), the first excited state is the symmetric superposition, and the second excited 
state is m z — |0). (c) Plots of the majority and minority spin projected density of states (PDOS) for CoH and CoH 2 . The 
difference in majority and minority spin spectral weights indicate that CoH has a total spin S — 1 and CoH 2 has S — 1/2. 
(d) Plot of the asymmetry between majority and minority PDOS for CoH adsorbed on N at the ri (left) and v (right) high 
symmetry points. 


els producing a non-negligible in-plane component of the 
anisotropy. 

To model the experimentally observed tunneling spec¬ 
tra and to determine the magnetic anisotropy we use a 
phenomenological spin Hamiltonian including the Zee- 
man energy and magnetic anisotropy: 

H = 9f x B B-S + DS 2 z +E(S 2 x -S 2 y ), (1) 

with g as the gyromagnetic factor, B the magnetic field, 
S = ( S x , S y , S z ) the total spin operator, and D and E as 
the axial and transverse magnetic anisotropy [3 H2 HI¬ 
ES]. Transport through the junction is calculated using a 
Kondo-like interaction <r • S between the tunneling elec¬ 
trons and the localized spin system, with a as the stan¬ 
dard Pauli matrices. We account for scattering up to 3 rd 
order in the matrix elements by considering additional 
exchange processes between the localized spin and sub¬ 
strate electrons of the form [27] (see SI): 

\jpo&- s. (2) 

To confirm the magnetic origin of the spectroscopic 
features in CoH and C 0 H 2 , we measure the differential 


conductance at magnetic fields up to 14 T normal to 
the surface. Figure (3 }d shows experimental spectra taken 
over one C 0 H 2 complex and Figure [3]: the model calcu¬ 
lations for the Kondo resonance. Applying an external 
magnetic field introduces Zeeman splitting to the spin- 
1/2 system (Figure [3 ^l) . At low magnetic fields, 2.5 T, 
the peak broadens and the differential conductance of 
the resonance is reduced. Increasing the field to 5 T, a 
clear splitting of the Kondo resonance is observed. For 
the highest fields, the degeneracy of the spin-1/2 state 
is effectively lifted, resulting in a strong reduction of the 
Kondo resonance and the appearance of an inelastic exci¬ 
tation gap. We can reproduce the peak and its splitting 
by our perturbative model (Figure [ 3 J 3 ) even though at 
high fields the peak-like conductance is weaker in the 
experimental data than expected from the model calcu¬ 
lation. This indicates that the Kondo temperature of the 
system lies close to the base temperature of our experi¬ 
ment. 

Increasing the external magnetic field has two effects 
on the spin-1 CoH; Zeeman splitting separates the steps 
and the ratio between inner and outer conductance step 
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Figure 3. Magnetic field behavior of C0H2 and CoH. (a) Left: Zeeman splitting of the spin-1/2 states of a C 0 H 2 
complex in magnetic field. Dots mark the energy differences as determined by least-square fits of the perturbation model to 
the experimental data in (b). The regression line corresponds to a gyromagnetic factor g — 2.0 ±0.1. Right: Sketch of the 
C 0 H 2 complex adsorbed on a N site, (b) Evolution of the differential conductance of a C 0 H 2 complex in an external magnetic 
field normal to the surface (B = 0, 2.5, 5,10 and 14 T; T = 1.4 K). (c) Simulated spectra using a 3 rd order perturbation model 
and a constant coupling to the substrate of —Jpo = 0.1 and g — 2.0. (d) Left: Sketch of the spin-1 CoH complex adsorbed on 
a N site. Right: State energy evolution in magnetic field along the out-of-plane anisotropy axis. Dots mark the experimentally 
determined step positions, full lines are the calculated eigenstate energies of the model Hamiltonian (see text) using magnetic 
anisotropy parameters of D — —4.8 meV, E — 0.6 meV, and g — 2.2. (e) Evolution of the differential conductance of a CoH 
system in an external magnetic field normal to the surface (B = 0,2.5, 5, 7.5,10 and 14 T; T = 1.4 K). (f) Simulated spectra 
using the parameter from (d) and — Jpo = 0.1. The 14 T spectrum is shown together with a 2 nd order perturbation theory 
model, i.e. —Jpo = 0 (dashed line), to highlight the necessity of 3 rd order contributions. Curves in (b, c) and (e, f) are shifted 
vertically for better visibility. 


height decreases (Figure^). At zero field, the ground 
and first excited states are a superposition of m z = |±1) 
and | —1) states, applying a magnetic field reduces the 
spin mixing and leads towards a |±1) ground and | —1) 
excited state. This accounts for the reduction of the in¬ 
ner step with increased magnetic field, as the transition 
between ground and first excited state becomes less prob¬ 
able because it would require a change in m z of two. Re¬ 
verting to a purely 2nd order simulation, large deviations 
are observed at both steps, evidence that coupling of the 
spin to the substrate conduction electron bath must be 
considered (Figure dashed line). The experimental 
data fits excellently when including 3rd order terms, i.e. 
assuming a finite — Jpo, an out-of-plane anisotropy axis, 
and g = 2.2 ± 0.2. 

Evaluation of more than 30 CoH shows no sharp dis¬ 


tribution of the anisotropy parameters D and E. A tran¬ 
sition of the main anisotropy axis into the surface plane 
occurs when 3 E > |D|, therefore we have only consid¬ 
ered complexes with a clear out-of-plane anisotropy de¬ 
termined by the criterion \D\ /3 E > 1.5; a representative 
spectrum with in-plane anisotropy is shown in the SI. 
By considering the values of — Jpo from our fits, we ob¬ 
serve a correlation between the magnetic anisotropy and 
coupling with the substrate, —Jpo- The red branch in 
Figure[4^ shows that as the substrate coupling increases, 
the axial magnetic anisotropy decreases. These results 
are in line with predictions that increased coupling shifts 
energy levels. The solid red line shows the best fit to 
our data and follows the trend D = Dq( 1 — a(Jpo) 2 ), 
where a is a constant describing the bandwidth of the 
Kondo exchange interaction. The shaded red region ac- 
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Figure 4. Influence of environmental coupling on CoH spectra, (a) Experimentally determined D and E (red and blue 
dots) parameters plotted versus the coupling strengths —Jpo. Full lines show the expected renormalization of D and E due to 
virtual coherences calculated with a Bloch-Redfield approach taking exchange scattering with the dissipative substrate electron 
bath into account. Shaded region shows the experimental uncertainty, (b) Computed differential conductance for different 
coupling strengths between the localized spin and the electrons of the substrate ranging from —Jpo = 0 to —Jpo = 0.2. At 
stronger couplings (—Jpo > 0) an increase of the outer step’s shoulder is expected concomitant with a reduction of the energy 
position of the outer step. This is equivalent to a reduced anisotropy energy D. (c) Schematic diagram showing the effect of 
exchange. When the exchange coupling, J, between the local spin and the conduction electron bath is weak, a large magnetic 
anisotropy, D, is observed (top). As exchange coupling to the substrate strengthens, the magnetic anisotropy is reduced driving 
the system closer to a Kondo state (bottom). 


counts for the possible range of a by considering an ef¬ 
fective bandwidth of ujq = 0.4 — 1.2 eV (see SI). For 
0.1 < — Jpo < 0.2 the variation in magnetic anisotropy 
fits exceptionally well, but for small values of —Jpo, some 
spread in the axial anisotropy is observed. These fluctua¬ 
tions are not accounted for in our model and indicate that 
for small — Jpo additional factors such as strain or defects 
may contribute to the magnetic anisotropy. While the 
axial anisotropy shows clear dispersion, the transverse 
anisotropy is essentially constant (Figure |4^i, blue). 

Figure^ shows the influence of — Jpo on the tunnel¬ 
ing spectra calculated using a Bloch-Redfield approach 
to incorporate virtual correlations between the ground 
and excited states due to the coupling with the dissipa¬ 
tive spin bath in the substrate assuming a flat density 
of states and an effective bandwidth of ujo = 1 eV (see 
si) ®mm- As —Jpo is increased, virtual correla¬ 
tions lead to renormalization and reduce the level split¬ 
ting. This is observed experimentally as a reduction of 
the axial magnetic anisotropy. Furthermore, higher order 
scattering processes in the tunneling influence the con¬ 
ductance leading to an enhanced shoulder at the outer 
energy step that changes the contours of the spectrum 
(see SI). The symmetric peaks shift towards zero bias 
as — Jpo increases indicating that correlations drive the 
anisotropic spin-1 system closer to the Kondo state. Fig¬ 
ure & schematically depicts the observed trend, when 
the spin is weakly coupled to the conduction electrons 
the magnetic anisotropy is stabilized. Increasing the ex¬ 
change interaction introduces correlations between the 
excited spin states and the conduction electrons, leading 


to a net reduction in the magnetic anisotropy. 

In conclusion, our results show that the Kondo 
exchange interaction modulates the magnetic anisotropy 
of single spin CoH complexes. The role of exchange was 
quantitatively determined by exploiting the corrugated 
h -BN moire structure. In conjunction with 3rd order 
perturbation theory simulations, we extracted the 
precise values of the spin coupling to the environment 
and its influence on the magnetic anisotropy. Kondo 
exchange must be considered an additional degree of 
freedom - beyond local symmetry, coordination number, 
and spin state - for spins connected to contacts. This 
parameter is non-local and therefore expected to be 
discernable at surfaces, in junctions, and perhaps in 
bulk SMM materials. 
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METHODS 

The Rh(lll) surface was prepared by multiple cycles 
of argon ion sputtering and annealing to 1100 K. On the 
final annealing cycle borazine (B 3 N 3 H 6 ) was introduced 
at a pressure of 1.2 x 10 “ 6 mbar for 2 minutes resulting 
in a monolayer h -BN film. Cobalt was deposited onto a 
cold, ^ 20 K, h -BN surface via an electron beam evapo¬ 
rator. 

Scanning tunneling experiments were performed on a 
home-built STM/AFM in ultra-high vacuum with a base 
temperature of 1.4 K and magnetic fields up to 14 T. 
All spectroscopic ( dl/dV ) measurements presented were 
obtained with an external lock-in amplifier and a mod¬ 
ulation voltage of 0.2 mV applied to the bias voltage at 
a frequency of 799 Hz. The tunneling setpoint before 
the feedback loop was disabled was V = —15 mV and 
I = 500 pA. For measurements on the same adatoms 
in different external magnetic fields the tip was retracted 
while the field was ramped and allowed to settle for max¬ 
imum stability. 

First principles calculations have been carried out in 
the framework of the density functional theory (DFT) 
as implemented in the VASP code |28j [29] . We use the 
projector augmented-wave technique m where the ex¬ 
change and correlation were treated with the gradient- 
corrected PBE functional as formalized by Perdew, Burke 
and Ernzerhof [31] . Hubbard U and J values were taken 
from self-consistent calculations and fitting to experi¬ 
ments to be U — J = 3 eV [32H34] . Full details are 
presented in the Supplementary Information. 
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